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*0 . A measurement of the integrated luminosity at the ep collider HERA is presented, ex- 

^ ! ploiting the elastic QED Compton process ep — > ejp. The electron and the photon are 

detected in the backward calorimeter of the HI experiment. The integrated luminosity of 
the data recorded in 2003 to 2007 is determined with a precision of 2.3%. The measurement 
is found to be compatible with the corresponding result obtained using the Bethe-Heitler 
l Vl ' process. 
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1 Introduction 



For particle collider experiments, the precise knowledge of the luminosity is essential for any 
type of cross section measurement. The instantaneous luminosity is defined as 

L = ^, (1) 

where / is the revolution frequency for the two colliding particles p\ and p 2 , n is the number 
of colliding bunches per revolution, and Ni (N 2 ) is the number of particles of type pi (p 2 ) per 
bunch. The effective cross section of the beams is A. The time-integrated luminosity C relates 
the cross section a PlP2 ^x of the reaction p\p 2 — > X to the number of events N PlP2 _+x expected 
in the time interval T by 

C= I Ldt= N p^ x . (2) 

Since it is difficult to monitor all beam parameters with a per cent level precision, in particular 
those defining the effective beam cross section A, the integrated luminosity is often determined 
by counting the number of observed events for a specific reaction pip 2 — > X with a well-known 
cross section. 

At HERA, the colliding beams are protons and electrons^. For the data taking period studied 
in this analysis, the proton beam energy is = 920 GeV and the electron beam energy is 
.Eg = 27.6 GeV. The reaction used to determine the integrated luminosity is the production of a 
radiative photon in elastic ep scattering, ep — > cyp. Depending on the phase space considered, 
this process is referred to as Bethe-Heitler (BH) scattering or QED Compton (QEDC) scattering. 
In the BH process JT), both the electron and the photon are emitted almost collinearly to the 
incident electron. The corresponding cross section is very large, 0(100 mb). Dedicated small 
angle detectors are used to record BH events. In contrast, for QEDC scattering [2], the particles 
have a sizable transverse momentum with respect to the incident electron and can be detected 
in the main detector. The momentum transfer squared at the proton vertex, t, is generally small. 
At very small momentum transfer \t\ -C 1 GeV 2 , elastic scattering dominates. At \t\ > 1 GeV 2 , 
inelastic processes are relevant and the reaction is sensitive to the proton structure. In addition, 
there are quasi-elastic contributions to the cross section, where the outgoing proton forms an 
excited state, like A + or N*, which then decays to a low mass hadronic system. Within the 
phase space considered in this analysis, the elastic QEDC cross section is O(50 pb). 

At HERA, the integrated luminosity is usually measured in the BH process, using dedicated 
detectors located at small angles. The advantage of this process is its very large cross section, 
thus negligible statistical uncertainties are achieved for small amounts of integrated luminosity. 
However, there are various sources of possibly large systematic uncertainty. For example, there 
may be inevitable acceptance limitations for the small angle detectors, caused by elements of 
the beam transport system which separates the BH photons and electrons from the circulating 
proton and electron beams. The acceptances of the photon and electron detectors may ex- 
hibit complex spatial structures and can vary in time as well. Another complication originates 
from synchrotron radiation emitted by the electron beam as it passes the focusing magnets sur- 
rounding the interaction region. Furthermore ep collisions can happen also outside the nominal 



1 In this paper the term "electron" is used genetically to refer to both electrons and positrons. 
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interaction region. These contribute to the BH measurement of the integrated luminosity, but 
must be corrected for when analysing cross sections with the HI main detector, which has a 
more limited acceptance as a function of the collision point position. 

In this paper, a determination of the integrated luminosity is presented, based on the elastic 
QEDC process, which is measured in the HI main detector. This method is insensitive to details 
of the beam optics. However, the smallness of the cross section leads to limited statistical 
precision, thus time-dependencies can not be resolved with high resolution by counting elastic 
QEDC events alone. Auxiliary measurements of other reactions with higher cross section may 
be used to monitor time-dependencies. 

Comparisons of BH and QEDC measurements at HERA have been performed previously 
J3]|. The inelastic QEDC process also has been measured at large \t\ flU. The data available for 
the elastic QEDC analysis described in the following were recorded with the HI detector in the 
years 2003 to 2007. 



2 HI Detector 

In the following, only those components of the HI detector are briefly introduced which are 
essential for the present analysis. A detailed description of the whole detector in its original 
configuration can be found elsewhere [5J. Components which were part of later upgrades are 
referred to here separately. The origin of the HI coordinate system is the nominal ep interaction 
point. The direction of the proton beam defines the positive z-axis (forward direction). Trans- 
verse momenta are measured in the xy plane. Polar (9) and azimuthal (ip) angles are measured 
with respect to this reference system. The pseudo-rapidity is defined as rj = — lntan(6'/2). A 
schematic view of the HI detector with signals from an elastic QEDC candidate event is shown 
in figure \T\ 

In the backward region —4.0 < r] < —1.4, a lead- scintillating fibre calorimeter [6] (SpaCal) 
is used for the identification and measurement of both the scattered electron and the scattered 
photon. The energy resolution for electromagnetic showers is a (E)/E ~ 7.1%/v/E/GeV© 
1% 0. The electromagnetic section of the SpaCal is read out in cells of size 4 x 4 cm in the 
xy plane, where the Moliere radius is 2.5 cm. The xy position of a shower is reconstructed as 
a weighted mean of the cell centres, the weights taken proportional to the logarithm of the cell 
energies [8j. After applying xy dependent corrections, the position resolution is about 3.5 mm 
for an electromagnetic shower in the energy range relevant to this analysis. 

The liquid argon (LAr) calorimeter covers the range — 1.5 < r\ < 3.4. Its electromagnetic 
(hadronic) section is equipped with absorbers made of lead (steel) plates. An energy resolution 
of a(E)/E ~ 11%/y/E /Ge V for electromagnetic showers and of a (E)/E ~ 50%/y/E/GeV 
for hadronic showers is obtained from test beam measurements. 

The central region of the detector is equipped with a set of tracking detectors (CTD). There 
are the two concentric central jet chambers (CJC), interleaved by a z chamber, and the central 
silicon tracker (CST) [9 |. The CTD measures the momenta of charged particles in the angular 
range 20° < 6 < 160°. The central inner proportional chambers (CIP) [fTOl are located between 
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the inner CJC and the CST. The five CIP chambers have a radial spacing of 9 mm, where the 
innermost layer is located at a radius of 15.7 cm. In if there is a 16-fold segmentation, whereas 
in z the segments have variable size, ranging from 1.8 cm in the innermost layer to 2.3 cm in 
the outermost layer. The CIP has an angular acceptance in the range 10° < 9 < 170°. In the 
backward region, the tracking is complemented by the backward proportional chamber (BPC), 
located directly in front of the SpaCal. 

The calorimeters and tracking detectors are located inside a large superconducting solenoid, 
providing a uniform field of 1.16 T strength. The return yoke of the solenoid is instrumented 
and serves as a muon detector. Upstream and downstream of the interaction point there are 
systems of scintillators (VETO), providing time-of-flight information. Timing signals from the 
VETO and the SpaCal were used during data taking to reject particles originating from non-ep 
interactions of the proton beam in the HERA tunnel. The luminosity system for measuring 
the Bethe-Heitler process consists of an electron tagger located at z = —5.4 m and a photon 
calorimeter located at z = —103 m. 



3 Signal and Background processes 

Monte Carlo event generators (MC) are used to predict event yields of signal and background 
processes. A GEANT3 iTTTTl simulation of the HI detector is performed for each generated event, 
where also the relevant time-dependencies such as changes to the detector setup and varying 
beam conditions are taken into account. Electromagnetic showers are simulated using a shower 
library [[8). After detector simulation, the events are passed through the same reconstruction 
algorithms as were used for the data. 

The QEDC signal is simulated using the COMPTON22 event generator [fT2l . This generator 
produces elastic, quasi-elastic and inelastic events. The elastic QEDC events are taken as signal, 
since their cross section only depends on QED theory and on the proton elastic form factors, 
thus having small uncertainties. Details are discussed in Section l54l The quasi-elastic events 
are treated as background and suppressed in the analysis, because their cross section depends on 
less precisely known parameters such as probabilities to produce excited nucleons. Similarly, 
the inelastic events are treated as background, because their cross section depends, for exam- 
ple, on parameterisations of the proton structure functions at very low momentum transfer. In 
COMPTON22, the fragmentation of quasi-elastic events is modelled using the SOPHIA pack- 
age |[T3l . whereas for inelastic events string fragmentation as implemented in PYTHIA QjJ is 
used. For the elastic QEDC signal, final state radiation from the electron has been included in 
the COMPTON22 event generator using the relevant PYTHIA routines. 

An important source of background is electron-positron pair production, ep — > epe~e + , 
simulated using the GRAPE event generator lfl~5Tl . Other background events originate mainly 
from various diffractive processes, namely deeply virtual Compton scattering (DVCS), diffrac- 
tive vector meson (VM) production and non-resonant diffraction. DVCS is modelled using 
the MILOU event generator |fl~6l. Diffractive VM production is simulated using the DIFFVM 
event generator [fTTl , where the production of p°, cu, <p, J/ip, ip' and T mesons is considered. 
For p° production, DIFFVM is modified such that decays to 7r°7 and 7/°7 are included. Non- 
resonant diffraction is simulated using the RAPGAP event generator [[T8~l . Background from 
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non-diffractive deep-inelastic scattering (DIS) is simulated using DJANGO |[T9ll and is found to 
be negligible. 



4 Event selection 

During data taking, events with electrons in the SpaCal were recorded if certain trigger criteria 
were fulfilled. The trigger required the presence of at least one localised energy deposit in the 
SpaCal calorimeter with energy E > 6 GeV. In addition, there were veto conditions to reject 
beam related background not originating from ep collisions. Data periods where the SpaCal 
trigger was not fully efficient are removed from the analysis, corresponding to about 6% of the 
total HI data sample recorded in 2003 to 2007. 

The elastic QEDC events are selected offline by requiring two clusters in the electromagnetic 
section of the SpaCal. A summary of the selection criteria is given in table [T] The transverse 
sizes of the SpaCal clusters are restricted to R\ og < 6 cm, where R\ og is calculated from the 
SpaCal cell centres using logarithmic energy weighting |JSj . The cluster energies are required 
to be larger than 2.2 GeV. In the range 30 < R < 72 cm of the radial distance from the beam, 
R, exactly two such clusters are required, whereas for 20 < R < 30 cm no cluster is allowed. 
The restriction in the number of clusters removes background and events with hard radiative 
photons. The restriction R > 30 cm on the two clusters ensures that the particles are within the 
CIP acceptance. 

Electron trajectories are reconstructed using the SpaCal cluster position together with posi- 
tion information from the CIP chambers. Hits in the CIP chambers are considered if they match 
the SpaCal cluster in (p. Adjacent hits are merged to CIP clusters in ip and z. This merging is 
done separately for each layer. A straight line fit of the CIP clusters and the SpaCal cluster in 
the rz plane is performed, where outliers are rejected. The coordinate r is the radial distance 
from the z axis, measured along the azimuthal direction given by the SpaCal cluster. After out- 
lier rejection, there are up to five accepted CIP clusters, corresponding to the five CIP layers. 
At least two accepted CIP clusters are required. Next, the centre-of-gravity of the CIP clusters 
in the rz plane is calculated. Finally, the CIP centre-of-gravity in rz , together with the SpaCal 
energy and the SpaCal position are used to reconstruct a helix trajectory in three dimensions, 
pointing back to the origin of the interaction. In the determination of this helix, the beam spot 
and beam tilt are also usecH. The direction of bending in the magnetic field is chosen assuming 
that the particle charge is equal to the charge of the beam lepton. The algorithm finally returns 
the origin of the interaction (CIP vertex) and the momentum vector at the CIP vertex. 

The electron and the photon are then identified, making use of the helix fit results. If there 
is no CIP vertex, the event is rejected. If there is only one SpaCal cluster linked to a CIP 
vertex, that cluster is taken as electron while the other cluster is taken as photon. The photon 
momentum vector is calculated from the photon cluster energy and a straight line trajectory 
pointing from the electron CIP vertex to the photon cluster position. If both SpaCal clusters 

2 The beam spot is defined as the average x and y position of interactions which take place at z = 0. The 
beam tilt is a slope correction for interactions at z ^ 0. These parameters were monitored regularly in short time 
intervals using the CTD. 



7 



Electromagnetic SpaCal clusters with transverse sizes R\ og < 6 cm and E > 2.2 GeV 
Exactly two clusters with radial distance to the beam axis 30 < R < 72 cm 
No additional cluster with 20 < R < 30 cm with energy E > 2.2 GeV 
At least one of the two clusters with CIP vertex 
Electron and photon identification 

Energies min(£' e , _E 7 ) > 7 GeV and max(E e , E 7 ) > 10 GeV 

Vertex position |;z vtx | < 35 cm 

Polar angles 155.9° < 9 e , # 7 < 169.5° 

Difference in azimuth 170° < Ay? < 190° 

Transverse momentum balance |P™ SS | < 0.3 GeV 

Energy in LAr at 9 < 10° is below 0.5 GeV 

No additional CTD track 

Table 1: Summary of selection criteria. 



are linked to CIP vertices, it is assumed that the photon has converted into an electron-positron 
pair while passing the material in front of the CIP detector. For that case, two hypotheses 
are checked in the reconstruction. First, one of the clusters is taken as the electron, and the 
photon momentum is calculated using the energy and position of the other cluster as described 
above. The difference in azimuth between the electron and photon candidate momenta, A</?i, 
is determined. Next, the particle hypotheses are interchanged and the corresponding difference 
in azimuth, Ay? 2 , is calculated. The particle assignment is done according to the hypothesis 
yielding an azimuthal opening angle closer to 180°. In the simulation, the mis-identification 
probability is 0.3% (16%) if one (both) SpaCal clusters are linked to CIP vertices. 

Once the electron and photon are identified, the z position of the electron CIP vertex, z vtx , 
is verified. Only events with \z vtx \ < 35cm are selected in the analysis. To further suppress 
background contributions the following cuts are applied: energy of the most (least) energetic 
particle greater than 10 (7) GeV, polar angles 9 ei 9 1 within 155.9° and 169.5°, difference in 
azimuth between 170° and 190°, modulus of the transverse component of the missing momen- 
tum smaller than 0.3 GeV. The missing momentum vector is calculated from the reconstructed 
electron and photon four-momenta. 

The inelastic background sources are further suppressed by using conditions on additional 
activity in the detector. Events are rejected if the energy in the forward part of the LAr calorime- 
ter, with polar angle 9 < 10°, exceeds 0.5 GeV, or if there are CTD tracks which can not be 
attributed to either the electron or the photon. The total number of elastic QEDC candidate 
events is 14277, after correcting for trigger efficiency effects. 

The efficiency of the reconstruction is determined for the following phase space defined for 
the generated elastic QEDC signal: the polar angular range of the generated electron and photon 
is restricted to 155.9° < 9 < 169.5°, the maximum fraction of incident electron energy carried 
away by initial state radiation is 0.35, and the momentum transfer squared at the proton vertex is 
limited to \t\ < 0.09 GeV 2 . For this generated phase space, the cross section is a gen = 54.8 pb. 
A total fraction of / vtx = 2.5% of the generated events have a simulated vertex position along 
the z coordinate, z gen , outside ±35 cm around the nominal vertex position, in order to be able 
to describe the observed longitudinal vertex distribution in the HI detector. Within the region 
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no |P™ SS | cut 


|P^ 1SS | < 0.3 GeV 


quasi-elastic QEDC 


6.84% 


2.93% 


inelastic QEDC 


7.02% 


1.51% 


elastic DVCS 


2.10% 


1.26% 


quasi-elastic DVCS 


0.55% 


0.16% 


ep — > ep e + e~ 


1.15% 


1.31% 


diffractive DIS 


2.78% 


0.53% 


non-diffractive DIS 


0.02% 


0.01% 


diffractive p° 


2.05% 


0.15% 


diffractive to 


0.43% 


0.03% 


diffractive 


0.29% 


0.02% 


diffractive J/ip 


0.20% 


0.05% 


diffractive tp' 


0.17% 


0.08% 


diffractive T 


0.02% 


0.01% 



Table 2: Background processes contributing to the elastic QEDC selection. The background 
fractions are given for both a selection with the |P™ SS | cut released and for the final selection. 



\Zgen\ < 35 cm the reconstruction efficiency is found to be e rec = 64.7%. The sources of 
inefficiency are investigated in the following. Losses of 12.4% originate from the cut in IP™ 1815 ! 
due to the limited experimental resolution. The requirement 30 < R < 72 cm on the radial 
SpaCal cluster position reduces the efficiency by 7.2%. This loss is related to a geometrical 
effect, such that particles originating from z gen > (z ee n < 0) and scattered at polar angles near 
155.9° (169.5°) are outside the allowed range in R. The track and calorimeter veto conditions 
contribute to the inefficiency by 3.4%, dominated by the restriction in forward LAr energy and 
the veto on additional SpaCal clusters. The LAr condition is not fully efficient due to electronic 
noise and overlap with non ep background. The inefficiency due to a third SpaCal cluster 
originates from events with hard final state radiation. Losses due to the other selection criteria 
are small (1.8% total). Finally, 10.5% are rejected by more than one selection criterion, where 
combinations involving the |P™ SS | condition dominate. 

A detailed breakdown of the different background sources contributing to the elastic QEDC 
selection as defined in table[T]is given in table |2] without and with applying the |P™ SS | cut. More 
than half of the background originates from inelastic and quasi-elastic QEDC processes. The 
|P™ SS | cut significantly reduces the background by about a factor of three. 

For measuring the integrated luminosity £°- EDC of collisions originating from the region 
\zgen\ < 35 cm, the following relation is used 

^ QEDC = (1 _ ^ jVeventCWbgr^ (3) 

where A^ even t is the number of QEDC candidate events observed in the detector, /b gr is the 
background fraction predicted by the MC simulation and a v i s = 36.4 pb is the visible QEDC 
cross section. The main contribution to cr vis originates from genuine QEDC production in the 
phase space region of this analysis and with an interaction vertex within the accepted region 
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(1 — /vtxVgen x e rec = 34.6 pb. Additional contributions from events outside the defined phase 
space or with an interaction vertex beyond the defined limits amount to 1.8 pb. 

Distributions of variables used in the selection procedure are shown in figure |2] Within 
uncertainties, the data are described by the prediction. Note that the prediction is normalised 
to the integrated luminosity £Q EDC as determined in the present analysis. The |P™ SS | distribu- 
tion, figure |2K is of particular interest as it shows a clear separation of background and signal. 
The analysis cut of 0.3 GeV is a compromise between inevitable systematic uncertainties due 
to the limited detector resolution, dominating at small |P™ 1SS |, and background contributions, 
increasing at large |P™ SS |. 

5 Systematic uncertainties 

Systematic uncertainties on the elastic QEDC measurement may be categorised as follows: 
experimental uncertainties, background uncertainties and QEDC theory uncertainties. The ex- 
perimental uncertainties originate from two sources, trigger and reconstruction efficiencies. A 
summary of the systematic uncertainties is given in table|3] The individual contributions are dis- 
cussed below. Additional time-dependent uncertainties are present in cases where the integrated 
luminosity determined in the present analysis is applied to subsets of the HI data. 

5.1 Trigger uncertainties 

The main trigger condition is based on calorimetric information in the SpaCal. It has an ef- 
ficiency of more than 95% for clusters with energies E > 6 GeV, rising to above 99.8% for 
energies E > 10 GeV. These efficiencies are verified using independent trigger conditions for 
a selection of DIS events with the scattered electron in the SpaCal. Both the electron and the 
photon from the elastic QEDC reaction create clear signals above the trigger condition thresh- 
olds, hence the trigger inefficiency on the QEDC selection is negligible. However, for certain 
time periods there were small regions opposite in (p with reduced trigger efficiency. This leads 
to an uncertainty of 0.02%. The other trigger conditions are related to timing signals from the 
VETO system and from the SpaCal calorimeter, designed to veto non ep background. These 
trigger conditions in conjunction with the varying HERA beam conditions cause inefficiencies 
of typically 1% for data taken up to the year 2005 and of typically 0.2% afterwards. These inef- 
ficiencies mainly originate from particles from beam related background recorded within gen- 
uine ep collision events. The veto inefficiencies are corrected for by applying time-dependent 
weights to the data events. The corresponding systematic uncertainty is 0.2%. 

5.2 Reconstruction uncertainties 

Reconstruction uncertainties originate mainly from the understanding of the SpaCal response 
to electrons and photons. The primary SpaCal energy calibration is done using electrons in DIS 
events ll20l . It corrects for time-dependent or spatial non-uniformities of the calorimeter re- 
sponse. However, the response of the SpaCal is slightly different for electrons and photons, 
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Experimental uncertainties 


trigger inefficiency 


0.2% 


SpaCal energy scale 


0.6% 


SpaCal energy resolution 


1.1% 


SpaCal position resolution 


0.3% 


CIP efficiency 


0.2% 


conversion probability 


0.3% 


alignment 


0.4% 


z-vertex distribution 


0.1% 


SpaCal cluster finder 


0.04% 


CTD efficiency 


0.03% 


LAr energy veto 


0.05% 


1.4% 


Background uncertainties 


non-elastic QEDC 


1.1% 


elastic DVCS 


0.3% 


quasi-elastic DVCS 


0.1% 


diffractive VM production 


0.1% 


non-resonant diffractive DIS 


0.2% 


ep — > epe + e~ 


0.1% 


1.2% 


QEDC theory uncertainties 


higher order corrections 


1.0% 


proton form factor (TPE parametrisation) 


0.1% 


proton form factor (experimental) 


0.2% 


size of generated signal sample 


0.3% 


1.1% 



Table 3: Systematic uncertainties on the determination of the integrated luminosity using elastic 
QEDC events. The different error sources are grouped into three categories: experimental, 
background related and theory related uncertainties. The error sources are described in detail in 
the text. 



mostly due to the presence of dead material in front of the calorimeter and due to final state 
radiation of the electrons. Furthermore, it is found that the primary calibration can be improved 
by correcting the energy response as a function of the transverse cluster size, i?i og . For the 
QEDC analysis, multiplicative calibration factors are applied to the SpaCal cluster energies for 
electrons, non-converted photons and converted photons, respectively. These factors are taken 
to have a linear dependence on i?i og . The corresponding parameters are determined by applying 
the double-angle calibration method to both the photon and the electron differentially in R\ og . 
Distributions of Pt/Pt,da, where Pt is the measured transverse momentum and Pt,da is the 
predicted transverse momentum, are investigated for the selection of QEDC events with the cut 
on the momentum balance, |P™ SS | < 0.3 GeV, relaxed. The predicted transverse momentum is 
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given by EM 

P T , DA = 2%( 1 -^^ + 1 -^lY\ (4) 
\ smfc'e sin 6^ / 

where 9 e and # 7 are the polar angles of the electron and the photon, respectively. For each 
Pt/Pt,da distribution, the position of the maximum is determined in a fit. The calibration 
parameters are finally determined from a linear fit as a function of R\ og . In data (MC), the 
energy response of the calorimeter to non-converted photons is found to be on average 3.5% 
(2.2%) higher than the response to electrons. For the event sample of converted photons, the 
energy response to photon candidates is 0.4% lower than the response to electrons, both for 
data and MC. In order to determine the systematic uncertainty, the energy scale of electrons and 
photons is varied separately by 0.5% each. The size of this variation covers possible systematic 
effects originating from the calibration procedure described above. Furthermore, a simultaneous 
variation of the electron and photon energy scale by another 0.5% is considered as systematic 
uncertainty, originating from the primary energy calibration ||20l . In total, all energy scale 
variations together contribute to the uncertainty on £Q EDC by 0.6%. 

In addition to the calibration factors, the energy resolution is determined from fits to the 
Pt/Pt,da distributions, however without dividing the sample into bins of i?i og . Figure [3] shows 
the distributions of Pr,e/ Pt,da and Pt^I Pt,da for electrons and photons, respectively. The 
distributions are peaked at 1, as expected after calibration. Near the peak, the distribution is 
more asymmetric towards smaller energies for electrons as compared to photons. This is at- 
tributed to final state radiation and energy losses in the material located in front of the calorime- 
ter. The original MC simulation (not shown in figure |3]) has deficits to describe both widths 
and tails towards lower transverse momenta. The effect exists for both electrons and photons. 
It is corrected for by applying an extra smearing of the reconstructed energies in the MC sim- 
ulation. An energy offset AE = (5 — t)E®/2 is subtracted, where 5 is a random number 
drawn from an exponential distribution, i.e. the probability density to find 5 > is given by 
f(5) = 1/r exp[— S/t]. By construction, AE ~ (5 — r) has an expectation value of (AE) = 0. 
This has the desired effect that the peak position of Pt/Pt,da is affected only little by the 
smearing. Two independent parameters r e and r 7 are foreseen to describe the expectation val- 
ues of the exponential probability distributions for electrons and photons, respectively. It turns 
out that both r e and r 7 take the same central value, r e = r 7 = 0.010. Figured shows that 
the distribution of Pr, e / Pt,da is described by the smeared simulation within a variation of the 
smearing parameter r e = 0.010 ± 0.005. Similarly, for photons, Pt^/Pt,da (Figure |3j>) is 
described within the variation r 7 = 0.010 ± 0.005. The r e and r 7 variations together cause an 
uncertainty on £Q EDC of 1.1%. 

The SpaCal cluster position resolution in MC is worse than in data. This effect has been 
identified using the difference in azimuth of the electron and the photon, Aip, shown in figure |4] 
The original simulation has a deficit at Aip near 180°, corresponding to a resolution worse than 
in data. In order to improve the description of data by MC, the reconstructed cluster positions 
in MC, x rec , are modified such that they are closer to the extrapolated SpaCal positions of the 
corresponding generated particles x gen . For the analysis, the positions x M c — (1 — f)x rec + fx seD 
arc used, where the constant is found to be / = 0.14. The uncertainty of / is taken as 0.05, 
resulting in an uncertainty on £Q EDC of 0.3%. The data are described by the prediction within 
that systematic variation, as demonstrated in figure @] 
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The CIP efficiency for electrons is determined in data and in the simulation using DIS events. 
It is found to be near 99% in data and near 99.5% for MC, varying as a function of the radial 
distance of the SpaCal cluster from the beam, R. A correction as a function of R is made by 
dropping a fraction of CIP vertices in the simulation. The CIP spatial resolution is adjusted 
using elastic p° production events in DIS. The CIP vertex, reconstructed from the scattered 
electron, is compared to the CTD vertex, reconstructed from the tx + tx~ pair. The conversion 
rate of photons in front of the CIP is underestimated in the simulation. In data, the conversion 
probability is around 32%, whereas the MC predicts 23%. This is corrected by mimicking 
conversion effects for a fraction of MC events with non-converted photons. For these events, 
extra CIP clusters near the expected position are added, and the energy response is scaled to 
match the expectation for converted photons. For estimating systematic effects, the three CIP 
related corrections described above are switched off one by one, and the resulting differences on 
£Qedc are ta k en as uncertainties. For the CIP efficiency correction the uncertainty is 0.2%. The 
CIP resolution tuning has negligible effect and the conversion probability leads to an uncertainty 
of 0.3%. 

The alignment of the SpaCal and CIP detectors is done using DIS events. The interaction 
vertex is reconstructed using tracks in the CTD, originating from the hadronic final state. Using 
hits in the BPC detector and the energy measured in the SpaCal, the electron trajectory is ex- 
trapolated to the CIP and SpaCal detectors. The alignment uncertainties are dominated by the 
uncertainty on the SpaCal z position. Systematic effects are estimated by varying the SpaCal z 
position by ±5 mm, resulting in an uncertainty of 0.4% on £ QEDC . 

The longitudinal vertex distribution is dominated by a Gaussian near z = with a width 
of approximately 10 cm, as can be seen in figure |2j. The longitudinal proton beam profile 
also exhibits prominent satellite peaks of similar width, leading to collisions in the HI detector 
near ±70 cm. In addition, there is an excess of collisions near 40 cm, as compared to a simple 
model which includes only collisions from the main bunch and from the satellites. For this 
analysis, the simulated vertex distribution is re-weighted such that the full interaction region 
is described. The difference of 0.1% in £ QEDC , obtained when using the simple beam profile 
model, is taken as a systematic uncertainty. The reconstructed z-vertex distribution after re- 
weighting is compared to the data in figure HJ. Good agreement is found. The regions of sizable 
systematic uncertainty due to the vertex re-weighting are visible. 

The identification of clusters in the SpaCal is checked by relaxing the R\ og < 6 cm condi- 
tion. Removing the R\ og condition results in a somewhat smaller number of selected QEDC 
events, because a third SpaCal cluster is accepted more often, leading to the rejection of the 
corresponding events. This procedure leads to a change in £Q EDC of 0.04% which is consid- 
ered as systematic uncertainty related to the cut in R\ og . The uncertainty on the CTD track 
reconstruction efficiency of typically 2% per track affects the analysis through the track veto, 
resulting in an uncertainty of 0.03% on £ QEDC . The veto on the energy in the forward part of 
the LAr calorimeter is checked by relaxing the veto condition to E < 1 GeV, resulting in an 
uncertainty of 0.05% on £ QEDC . 

5.3 Background uncertainties 

The normalisation of quasi-elastic and inelastic QEDC events predicted by the COMPTON22 
event generator depends mainly on parameters related to the nucleon excitation and on the 
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proton structure function parametrisation at low momentum transfer, respectively. These pa- 
rameters are not known very well. For this reason, the normalisation of the sum of these con- 
tributions, referred to as "non-elastic QEDC", is verified by investigating the vector sum of 
the electron and photon transverse momenta, P™ m = — P™ ss . The vector P™ m is decomposed 
into components parallel to (Pj!) and perpendicular to (P^r ) the electron transverse momentum. 
The distributions of Pj. and P^r are shown in figure |5] inside the analysis phase space as well 
as for |P™ 1SS | > 0.3 GeV with all other selection criteria applied. Both the parallel and the 
perpendicular components are described well. Outside the nominal analysis phase space the 
non-elastic QEDC contributions dominate at large p| or at large P?- The p| distribution is 
somewhat asymmetric for |P™ SS | > 0.3 GeV, because in contrast to photons the SpaCal re- 
sponse to electrons has tails towards low energies, as discussed in the previous section. The 
normalisation of the non-elastic QEDC contribution is tested by performing fits to either p| or 
Py for |P™ SS | > 0.3 GeV. The normalisation factors observed in these fits are compatible with 
the COMPTON22 prediction within 25%, which is used as normalisation uncertainty for the 
non-elastic QEDC processes. 

The DVCS cross section predictions obtained with the MILOU program are in agreement 
with recent HI measurements ||23l . Uncertainties of 20% for the elastic DVCS process and 50% 
for proton dissociative DVCS are considered. 

The elastic VM production rates are normalised using dedicated selections as close as pos- 
sible to the QEDC analysis. However, instead of requiring a photon in the SpaCal, a vector 
meson is reconstructed. The decays p° — > ir + 7i~, — > K + K~, J/tp — > £ + £~ (£ = e, p) and 
ip' — > £ + £~ are reconstructed from two oppositely charged tracks, detected in the CTD. The 
decay u — > 7r + 7r~7r° is reconstructed from two charged tracks and one or two neutral calorime- 
ter clusters. The decay T — > e + e~ is reconstructed using a sample of photoproduction events, 
where an e + e~ pair from the T decay is reconstructed in the SpaCal, one of the SpaCal clusters 
matched with a central track, and the scattered electron is outside the acceptance of the HI de- 
tector. The following normalisation uncertainties are found: 20% on p° and <f) production, 50% 
on J/ip, and 100% on u, ip' and T. Possible contributions from p(1450) production are covered 
by the p° normalisation uncertainty. 

The rate of non-resonant diffractive events, simulated using RAPGAP, is normalised using 
a selection of low multiplicity final states, where the electron is reconstructed in the SpaCal and 
one up to three additional charged or neutral particles are found. The uncertainty is estimated 
to be 30%. 

For the QED processes modelled by GRAPE, an uncertainty of 10% is assumed, taking into 
account possible higher order effects. 

The uncertainties on the various background samples originating from the finite sizes of the 
generated event samples are negligible as compared to the uncertainties discussed above. 

5.4 QEDC theory uncertainties 

Uncertainties to the elastic QEDC cross section arise mainly from two sources: higher order 
corrections and the knowledge of the proton form factors. 
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In the original COMPTON22 event generator, higher orders are simulated in the peaking ap- 
proximation [24] . Improved higher order corrections have been calculated [1231 using a photon 
radiator [|26l . For the purpose of this analysis, the COMPTON22 events are assigned weights, 
determined such that the cross section predicted by the photon radiator method is reproduced. 
The difference of 1.0% to the original COMPTON22 prediction is taken as systematic uncer- 
tainty due to higher order effects. The elastic QEDC cross section also depends on the proton 
electric and magnetic form factors. In the original COMPTON22 generator, only a simple 
parametrisation of the form factors, using one parameter, is implemented. For this analysis, 
recent form-factor fits of elastic ep scattering data ll27l are taken into account, using an event 
weighting technique. The form-factor parameterisations [|2~7l are corrected for the effects of 
two-photon exchange (TPE), but parameterisations not including TPE corrections are also pro- 
vided. In the COMPTON22 generator, TPE contributions are not included when calculating the 
cross-section. For this reason the form factor parameterisations not including TPE corrections 
are used to calculate elastic QEDC cross sections [|28l . When computing the COMPTON22 
cross section with the TPE corrected form factors, the analysis result changes by 0.1%. This 
difference is included as systematic uncertainty. Experimental uncertainties on the form fac- 
tor parametrisation are also considered. Such uncertainties are available with the parametrisa- 
tion [29 J which includes TPE corrections. It has been verified that the difference between ||2~9~1 
and ll2~7ll is completely negligible for the purpose of this analysis if the TPE corrections are in- 
cluded. The experimental uncertainties on the form factors Ge(|£|) and Gm(|£|) approach zero 
for \t\ — > 0, because the parameterisations enforce Ge(0) = 1 andGM(O) = \x p . At \t\ = 0.007, 
which is the average momentum transfer predicted for the events selected in this analysis, the 
uncertainties on Ge (Gm) are 0.1% (0.2%). The elastic QEDC cross section at fixed \t\ is given 
by a linear combination of Cr#(|i|) and G 2 M (\t\), where the G\ contribution dominates at small 
\t\. When propagating the |t| dependent parametrisation uncertainties on G £ and Gm to the lu- 
minosity measurement, an uncertainty of 0.2% is found. In addition to the theory uncertainties 
related to higher order corrections and proton form factors , discussed above, there is a statistical 
uncertainty on the theory prediction originating from the finite size of the generated signal event 
sample. It amounts to 0.3%. 

In figure [6] the distribution of the variable (E — p z ) / (2E®) is studied. This variable is 
calculated from the sum of the four-momenta of the electron and the photon. The distribution 
of this variable is expected to peak at 1. The tail to small (E — p z )/(2E®) originates from 
initial state radiation, whereas values larger than 1 show up due to resolution effects. The data 
are described within the systematic and statistical uncertainties for both small and large values 
of (E — p z ) I (2£g). As expected, the peak region is dominated by experimental uncertainties, 
whereas the region of small (E — p z ) / (2E®) is dominated by uncertainties of the QEDC cross 
section. 



5.5 Time-dependent uncertainties 

In order to apply the integrated luminosity £Q EDC to other physics analyses, possibly using time 
restricted HI data sets, a luminosity calculation differential in time is required. This is achieved 
using DIS events measured in the SpaCal. The DIS selection follows the selection described 
in 11201 but is restricted in electron polar angle to the range 167° < 9 e < 172° such that the 
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expected event yield is most insensitive to changing beam conditions, in particular to the average 
position of the interaction vertex in z. In addition, the electron energy range is restricted to 15 < 
E < 25 GeV and the electron transverse cluster size to i?i og < 4.5 cm. The longitudinal vertex 
position, measured in the CTD, is restricted to be within ±35 cm around the nominal interaction 
point. The DIS event counts for each rurj^] are used to define relative integrated luminosities 
of the runs, and the overall normalisation is taken from the QEDC analysis. The statistical 
uncertainty of the DIS selection is negligible, but it has a time-dependent systematic error of 
1.5%. This uncertainty originates mainly from the SpaCal trigger and vertex reconstruction 
efficiencies ll2~0l . Figure [7J shows the results of the elastic QEDC analysis performed in bins of 
about 25 pb _1 , normalised to the global QEDC analysis with the DIS yield corrections applied. 
Four data taking periods, corresponding to distinct configurations of the HERA machine or the 
HI detector, are indicated. The HERA machine has been operated with e + p beams for periods 
(I) and (IV) and with e~p beams for period (II) and (III). The HI detector has been opened 
for the repair or upgrade of various components between period (I) and (II) as well as between 
period (II) and (III). The two methods of measuring differential in time are in good agreement, 
taking into account the statistical fluctuations of the time-dependent QEDC analysis and the 
time-dependent systematic uncertainties of the DIS yield method. 



6 Results 

The integrated luminosity of the data collected in the years 2003 to 2007 is determined using 
elastic QED Compton events. For the data sample as used in this paper, an integrated lumi- 
nosity of £Q EDC = 351.6 ± 8.0 pb -1 is measured. The statistical uncertainty amounts to 0.8%, 
whereas the total systematic error is 2.1%. The integrated luminosity is in agreement with the 
Bethe-Heitler measurement, £ BH = 338.9 ± 10.2 pb~ x . The corrections needed to measure the 
integrated luminosity of arbitrary time restricted data samples, such as samples comprising only 
e + p or only e~p beams, induce a further uncertainty of 1.5%. 



7 Conclusions 

The elastic QED Compton process is used to determine the integrated luminosity of the HI data 
taken in the years 2003 to 2007. The systematic uncertainties are about equally shared between 
experimental uncertainties, understanding of the elastic QEDC cross section and understand- 
ing of the background to the measurement. The statistical uncertainty is small compared to 
the systematic uncertainties. The new measurement method presented in this paper allow a 
determination of the integrated luminosity with a precision of 2.3%. 

3 H1 data is grouped into runs, where new runs are started whenever data taking conditions changed. A run 
typically spans about 30 minutes of data, equivalent to an integrated luminosity of about 30 nb _1 . 
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Figure 1: Elastic QEDC candidate event observed in the HI detector. The HI detector com- 
ponents most relevant to this analysis are indicated. The approximate electron and photon 
candidate trajectories are shown. 
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Figure 2: Distributions of variables used to select elastic QEDC events: (a-e) kinematic quanti- 
ties of the selected electron-photon pair and (f) the z coordinate of the position of the interaction. 
The kinematic quantities are (a) the minimum polar angle, (b) the maximum polar angle, (c) the 
minimum energy, (d) the maximum energy and (e) the modulus of their total transverse mo- 
mentum. The data are shown as black dots with the statistical uncertainties indicated as vertical 
bars. The simulation including background, normalised to the integrated luminosity determined 
in this analysis, is indicated as a solid line, with the systematic uncertainties attached as shaded 
area. Also shown is the contribution from background. The hatched areas are excluded by the 
selection criteria. 
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Figure 3: Distributions of the ratio of measured to predicted transverse momentum for (a) elec- 
trons and (b) photons. The predicted transverse momentum Pt,da is calculated using the double 
angle method. The data are shown as black dots. The simulation including background, nor- 
malised to the integrated luminosity as determined in this analysis, is shown as a solid line, 
with the systematic uncertainty originating from the limited knowledge of the energy resolution 
attached as a shaded area. Also shown is the contribution from background processes. 
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Figure 4: Distributions of the difference in azimuth between the electron and the photon. In 
(a) the event counts are shown, whereas in (b) the ratio to the expectation is drawn. The data 
are shown as black dots with statistical uncertainties indicated as vertical bars. The simulation 
including background, normalised to the integrated luminosity as determined in this analysis, is 
shown as a solid line with the systematic uncertainty originating from the limited knowledge of 
the position resolution attached as a shaded area. The distribution predicted by the simulation 
prior to adjusting the position resolution is shown by the dashed line. 
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Figure 5: Distributions of components of the photon plus electron transverse momentum sum, 
P™ m : (a) and (c) the component perpendicular to the electron transverse momentum, (b) and (d) 
the component parallel to the electron transverse momentum. The upper row, (a) and (b) shows 
the distributions inside the analysis phase space, the lower row (c) and (d) shows the distribu- 
tions for |P™ SS | > 0.3 GeV. The data are shown as black dots with the statistical uncertainties 
indicated as vertical bars. The simulation including background, normalised to the integrated 
luminosity as determined in this analysis, is indicated as a solid line, with the systematic uncer- 
tainties attached as shaded area. Also shown are the contributions from non-elastic QEDC and 
from other background sources. 
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Figure 6: Distribution of the variable (E — p z ) /2Eq calculated from the sum of the electron 
and photon four-momenta. In (a) the event counts are shown, whereas in (b) the ratio of data 
to expectation is drawn. The data are shown as black dots with the statistical uncertainties 
indicated as vertical bars. The simulation including background, normalised to the integrated 
luminosity as determined in this analysis, is indicated as a solid line, with various components 
of the systematic uncertainties attached as shaded areas. 
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Figure 7: Integrated luminosity measured from elastic QEDC events in bins of approximately 
25 pb -1 , divided by the integrated luminosity derived from the QEDC analysis on the full sam- 
ple with time-dependent corrections applied. The statistical uncertainties of the binned QEDC 
analysis as well as the uncertainties of the time-dependent corrections, here applied to four data 
taking periods (I)-(IV), are indicated. 
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